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Abstract 



We study logarithms of the form ln(m 9 /m&) which arise in the inclusive 
semileptonic decay of a bottom quark to a quark of mass m„. We use the renor- 
malization group to resum the leading radiative corrections to these terms, 
of the form m^a™ ln n (m q /mb), m^a™ +1 ln n (m q /mb) and m q a™ ln n+1 (m q /mb). 
The first two resummations are trivial, while the latter involves a non-trivial 
mixing of four-fermi operators in the l/m^ expansion. We illustrate this tech- 
nique in a toy model in which the semileptonic decay is mediated by a vector 
interaction, before treating the more complicated case of left-handed decay. 
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I. INTRODUCTION 



The inclusive semileptonic decay rate of a hadron containing a single bottom quark may 
be written as a power series in Aqcd/^6- The leading term in this expansion is simply the 
width for the underlying quark- level process b — > qlv + qglv + . . ., which may be expanded 
in powers of a s (mb), 



G 2 F \ V q b\ 2 
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r (o)/. \ . as(m b ) {1) fa s (m b )\ (2) 



The coefficients r( n ) are functions of the scaled final quark mass, m q = m q jm b . When 
the process is computed at the parton level, the masses arise in limits of the phase space 
integration, and hence the masses which appear are the perturbatively defined pole masses, 
m q = mP° le . Taking the leptons to be massless, the tree level term is 

T (0) (m q ) = l-8mj-24mjlnm 9 + 8mj-m^. (1.2) 

The full expression for T^(m q ) has been computed analytically fl|], and is quite lengthy; 
expanding the result in powers of m q gives 

„ m , A x 25 2tt 2 /136 „„, A \ . 2 64tt 2 a „ 
Y w (m n ) = : f — + 321nm mi H mi 
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^\nrh q + ^ In 2 m q \m* + 0(mf). (1.5) 
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At order a 2 , only the graphs containing gluon vacuum polarization have been calculated, 
and only numerically. For example |2[ 

T (2) (0) = (-3.44/5 + Cl ) , 
r (2) (0.37) = (-0.75/? + c 2 ) , (1.6) 

where C\ and c 2 denote terms not proportional to (3q = 11 — 2nf/3. 

Since inclusive semileptonic bottom decays provide a means of measuring the CKM 
mixing angle \V c b\, it is useful to have as much information as possible about the size of the 
higher order corrections to T. In this paper, we use the operator product expansion and the 
renormalization group to sum to all orders leading logarithms of the form m^o;"ln n+1 m q , 
for n > 0, as well as terms of the form m 2 a"ln n m ? and m^a" +1 ln n m g . We will show that 
these corrections are straightforward to calculate using the renormalization group. 

We will see that the resummation of the "phase space" logarithms m q a™ ln n+1 rh q is 
particularly interesting, and it is to them that we will pay the most attention in what 
follows. However, because of the prefactor rht, these terms are not dominant as m q — ► 0, 
or in any other limit of the theory. In fact, they are smaller, in principle, than uncomputed 
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terms of the form m^a™ -1 log n rh q , since m q \ogrh q — > as rh q — > 0. On the other hand, since 
m c is not particularly small, these terms may be numerically significant for b — > c decays. 

Unfortunately, by the same token rh c ~ 0.37 is such a poor expansion parameter^ that 
the terms which we can compute using the renormalization group do not dominate those 
which we cannot compute as easily, and so these results may not be used directly to estimate 
the size of the higher order corrections to b — > X c iv decays. This is clear from examining 
the sizes of the various terms which contribute to 1^°) and F^ 1 ': 

r (0) (0.37) = 1 - 1.095 + 0.448 + 0.021 - 0.0004 + . . . = 0.372 (1.7) 
r (1) (0.37) = -2.41 - 1.85 + 10.67 - 8.06 + 1.46 - 0.372 - 0.005 + . . . = -0.568. 



where the order of the terms is the same as in Eqs. (|1.2|) and (|1.3| ), and we have included terms 
up to (D(mj?) in the expression for Significant cancelations occur in both expressions 
between terms of different order in m c ; in particular, there is a large cancelation between 
the 0{a s m i c ) and 0(a s m^,) terms. Similarly, we will find when expanding our resummed 
results that there is a large contribution (larger than the tree level rate!) at O (o^m^ In m c ). 
In analogy with the one-loop expression, we might expect a large cancelation between this 
term and the order a 2 m^lnm c term. However, this latter term is down by two powers 
of lnm c relative to the terms which we are resumming (requiring a three-loop anomalous 
dimension to resum), and so we have not calculated it. The motivation for our analysis is 
the insight it will afford us into the origin of a variety of higher order terms in the expression 
for the semileptonic width, rather than in any reliable estimate of the true size of higher 
order corrections. 

We will use the operator product expansion (OPE) and the heavy quark effective theory 
(HQET) in our analysis. The application of OPE techniques to inclusive semileptonic 
heavy quark decays was suggested in Refs. P,[10|, in which two distinct, but ultimately 
equivalent, methods were introduced. The two approaches differ in the treatment of the 
leptons in the final state. Since the leptons interact only weakly and electromagnetically 
with the quark currents which mediate the hadronic decay, there is freedom to integrate 
over their momenta at various stages of the calculation. 

Let us consider the decay B(P B ) — ► X C (P X ) + (£v)(q). This process is mediated by a 
term in the weak Hamiltonian, 

n w = ... + cYi - l5 )bh,(i - l5 )u • (1.8) 

In the approach of Ref. ||, the Hamiltonian is explicitly factorized into a product of a quark 
current, and a lepton current Jf. Then the differential rate dr/dg 2 dg • v is written as 



dr 



u - x -' v ~ L lxu (q 2 ,q-v)W^(q 2 ,q-v), (1.9) 



dq 2 dq ■ v 



1 For our numerical results, we use the HQET relation rh c = m£)/mB + 0(l/mi, c ), where Ttlm = 
(tom + 3m* M )/4 is the spin-averaged meson mass. This gives m c = 0.37 instead of the more 
commonly used value of 0.3. 
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FIG. 1. Typical diagrams contributing to ImTjJ^, J u }- 



where = Pj^/rriB is the velocity of the B meson and L^ v is the spin summed lepton 
tensor (L^ oc (q^qu — g^q 2 ) for massless leptons.) The nonperturbative hadronic tensor 
W^ v is related via the optical theorem to the imaginary part of the forward scattering 
amplitude [|9|,|l0fl, 

= J2( B \ J h l^)(X c | r h \B){2-nf5\P B - P x - q) 

x c 

= -2lm(B\i [ dxe^T^ix), J£(0)} \B) . (1.10) 

The time-ordered product is then written via an operator product expansion as a power se- 
ries in a s (m b ) and l/m b , as illustrated in Fig. [|. The resulting expression for the differential 
rate consists of a series of delta functions and derivatives of delta functions at the thresh- 
old for c quark production, followed by a cut in the complex q ■ v plane corresponding to 
gluon bremsstrahlung. When integrated over the appropriate variables, this yields a sensible 
prediction for differential decay widths as well as for the total semileptonic width. In this 
approach, the factor of m\ In m c in Eq. (|1.2| ) arises from the integration over the phase space 
variables q and q ■ v. It is important to note that this term is not related to the running 
of the operators in the OPE between /i = m b and fi = m c , since this running is performed 
before any phase space integrals are performed. 

By contrast, in the approach of Ref. Jl0|] the OPE is performed on the expression for the 



total, rather than the differential, rate. Once again, the width is written as the imaginary 
part of the forward scattering amplitude, 

T B ^x c eu ~ E 1mA(B -> X c ev -> B) 
x c 

~\m(B\T{H ] w ,H w }\B). (1.11) 

This version of the time ordered product is illustrated in Fig. ^ While this approach 
is completely equivalent to the other, performing the OPE after the integration over the 
lepton momenta gives us more insight into the origin of the terms involving lnm c . 

When the OPE is performed at the renormalization scale /x = mj, the time ordered 
product in Eq. ( |1 . 1 1| ) is written as a sum of local operators, 

\mT{T0 w ,Hw} — ► ai{m b )bb + — {a 2 i{m b )rh 2 c bb + a 2 2(m b )b (iD) 2 b 

+ a 2 ^m b )ba^G^ b] + O (l/rnf) (1.12) 
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FIG. 2. Diagrams contributing to lm.T{H w ,Tiw}- 
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FIG. 3. Contribution to the coefficient function of bTiccT2b in the OPE of lm.T{H w , T~tw}- 

(note that only the imaginary part of the time ordered product is needed). At the matching 
scale there are no factors of lnm c in the coefficient functions a^(m&); these logarithms are 
infrared effects which are contained in the matrix elements of local operators. However, 
since none of the operators at order l/m| contains explicit c quarks, none of their matrix 
elements depend on m c at leading order in a s ; therefore there is no term proportional to 
m 2 c In rh c in the tree level expression ( |1.2|) for the semileptonic decay rate. The first operators 
of interest which contain explicit c quarks are four quark operators of the form bTiCcT2b. 
These arise in the matching due to the graph in Fig. |^ and are of relative order 1/mf. 

In the next section we will consider a toy model in which the b quark decays via a vector 
current. In this model, the operator hhcc arises at order l/m\ in the OPE; its matrix 
element between b quarks is given by the diagram in Fig. |], giving a contribution of order 
m 3 . In m c to the total inclusive rate. For the physically relevant case of a left handed current 
the corresponding operator is hh c^c, for which the graph in Fig. § vanishes. This explains 
the lack of a term of order m|! lnm c in the total inclusive rate ( |l-2j ). In this case, the first 
logarithm of rh c arises at order 1 / , in the matrix element of the operator rhjihcc, giving 
the term of order m^\nrh c in the inclusive decay rate. 

Rather than leave these logarithms in the matrix elements of local operators, it is con- 
venient to scale the theory down to the renormalization point fi = m c , at which point the c 
quark is integrated out of the theory. Below this scale matrix elements can no longer depend 
on m c ; all such dependence has been transferred to the coefficient functions in the OPE. By 
including the leading QCD corrections in the renormalization group equation, the complete 
series of leading logarithms of the form m™a™ ln n+1 m c may be resummed, 



where n = 3 for a vector current and n = 4 for a left handed current. This calculation 
is technically more complicated for a left-handed current than a vector current, since it 



c m t 



" In rh c — > m\ 



c In rh c + C\a s In 2 fn c + c 2 oi 2 s In 3 rh c + . . . 
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FIG. 4. Contracting the charm quark fields in a four fermion operator gives a ln(m c /m&) 
contribution to the semileptonic width. 

involves the renormalization of the complete set of dimension seven operators. Therefore we 
will warm up in the next section with the simpler case of a vector current, before proceeding 
on to the realistic decay. 



II. DECAYS VIA A VECTOR CURRENT 

We consider the decay mediated by the hadronic current cry^b, coupled to the usual 
left-handed massless leptons. At tree level, the decay width of the b quark is given by 

IV = rVo[l - 2m c - 8m 2 c - 18m 3 c + 18m 5 c + 8m 6 c + 2m\ - m 8 c 

— 24ml 1°S — 24m^ log m c — 24m\ log m c ] , (2.1) 

where Ty = ml/192n 3 . The total decay rate may be written via the optical theorem in 
terms of the imaginary part of the forward scattering amplitude. Integrating explicitly over 
the leptons, this may be recast as the expectation value of the time-ordered product of the 
hadronic current and its conjugate, 

Tv= 2^I dq6 ~ l9 X {Bl T ih"c(x),cYK0)} \B) x i- (q^ - q 2 g^) . (2.2) 

The integral is taken over physical values of the total lepton four-momentum q^. 

We now develop the time-ordered product in an operator product expansion. The mo- 
mentum transfer is of order over almost the entire region of integration, so we organize 
the expansion in inverse powers of rrib rather than in inverse powers of q 2 . Simultaneously, 
we will expand the ordinary b quark field in terms of the mass-independent HQET field h, 
defined by [|J 

h(x) = — ^— exp(im b v ■ x) b(x) . (2.3) 

Here t> M = p^/rrib is the four- velocity of the b quark, which is fixed in the limit — » oo [Q. 
The Dirac matrix (1 + f)j2 projects onto the heavy quark part of the field operator, so h is 
a two-component, rather than a four-component, object. The exponential factor cancels out 
the large "on-shell" part of the b quark momentum. This procedure will make all dependence 
on the heavy quark mass mj, explicit in the operator product expansion. 
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We will order operators in inverse powers of m&. The operators of dimension less than 
six can have no more than two fermion fields, which must be h and h, since we are interested 
in the heavy quark decay process. The c quark fields are contracted; the integral over the c 
quark momentum p^ is equivalent to the integral over the lepton momentum = p^ 1 — p^ 
in Eq. (|2.2|). The leading operator, then, is of dimension three: 



ml — 



hh = T vo hh. (2.4) 
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The matrix element of hh may be expanded in inverse powers of [0 , 

(B\hh\B) = 2m B [l + 0(l/ml)] , (2.5) 

where all 1/m" corrections are independent of m c . Hence the leading operator reproduces 
the m c = decay rate IVo- 

The next operator in the expansion in 1/my, is of dimension four, 

— 2Yvo m chh. (2.6) 

The factor of m c is an explicit part of the operator, including its dependence on //. The 
operator hh is a conserved current in the HQET and does not run |TTJ . Since the operator 



product expansion is performed at the scale fi = m&, it is m c (mb) which appears in the 
matching. Since mb does not run below fi = nib, neither does the combination m c (nib)h h. 
However, leading logarithms are generated if the decay rate is expressed, as it usually is, 
in terms of m c (m c ) or m^ ole . All of these logarithms are resummed if we leave m c (m&) 
unexpanded. The situation is analogous for the dimension five operator 

-8T V0 m 2 c hh. (2.7) 

Matching at \i — nib yields m c {vrib) in Eqs. ( |2.6| ) and ( |2.7| ), which resums the leading loga- 
rithms of the form a n s log" m c in these terms. We note that there are other operators which 
arise at dimension five and higher, such as h(iD) 2 h and ha^G^h, but the matrix elements 
of these operators are proportional to QCD scale quantities such as Ai and A2 and do not 
yield powers of 1 12fl . They are included in an analysis which treats the nonperturbative 



power corrections to the parton model decay ||13 HT7jl . 

At dimension six we first encounter the four-fermion operators of the form hTiccr 2 h, 
which give rise at tree level to the term proportional to ml logm c in the total decay rate. We 
will now compute the leading logarithmic improvement of these terms. Keeping the leading 
powers of mb in the lepton tensor (g M g„ — q 2 'g / / 37r ■, the operator product expansion yields 
the dimension six operator 

2 

h^cc^h x -^(v^ - fa) ■ (2-8) 

With the identity if)h = h on the HQET field, this reduces to 

2 

JYI r 

— - hcch - h^ccjah . (2.9) 

37T L - 1 
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These operators are matched onto at the scale /x = m&. Closing the charm quark loop, 
they mix with the operator w? c hh at order a°, yielding the "phase space" logarithm 
— 24rVo ? ^c log^c- To resum the leading logarithms, we must consider the QCD renor- 
malization of the dimension six operators. However, we do not need the QCD correction 
to the mixing of the four quark operators with the quark bilinears, which is subleading and 
will only produce terms in the rate proportional to rn? c a n s log n m c . 

The renormalization of the four quark operators simplifies considerably if we first apply 
a Fierz transformation to bring them into the form h£\h~cT<iC We use the 577(3) identity 

6u6 kj = 2Tt j T& + ±6 ij 6 kl , (2.10) 

on the color indices, and include a factor of —1 from the exchange of the fermion fields h 
and c. The Fierz transformation yields a linear combination of operators of the form 

Ol = hT a hcT a c and O s a = hT a T a hcT a T a c , (2.11) 

where T a = 1, 7 5 , 7 M , 7^7 5 , a^. However, because the heavy field h has only two components, 
there are actually just four independent scalars of the form hTh. It is straightforward to 
show, then, that the allowed Dirac structures for the dimension six operators are 

hhcc, 
hhcrfc , 

h^ 5 hc^ 5 c, (2.12) 
7i7 M 7 5 /ic^7 M 7 5 c. 

With both singlet and octet color structures, we find a set of eight operators. 

In the HQET, the coupling of a gluon to a heavy quark line is given by the Feynman rule 
v ■ A, where is the gluon field [|J. In v ■ A = gauge, in which the gluon does not couple 
directly to the heavy quark, we see that the renormalization of 0\ and O® is restricted to 
the cTc part of the operator. Since for m c = the addition of a gluon loop to cTc generates 
an even number of Dirac matrices and no 7 5 , none of the eight dimension six operators 
mix with each other under renormalization. Hence the running of these eight operators is 
multiplicative, except for the color octet hT a hcf)T a c, which mixes via "penguin" diagrams 
with the flavor SU(3) singlet operator hT a h'q i fT a q,^ with % summed over u,d,s. 

The dimension six operators mix with quark bilinears at order a° s by contracting the 
charm quark fields as in Fig. However, only the color singlet, scalar-scalar operator 

6 = hhcc (2.13) 

has a color and Dirac structure such that this mixing, to rh^h h, is nonvanishing. Hence it is 
the only operator which we must consider. The Fierz transformation of the operators ( |2.9|) 
yields the coefficient function 

Ce(m b ) = ^f. (2.14) 
Cg(/i) runs according to the renormalization group equation 



S 



fi-^-CM = j 6 C 6 (ii) , 
where a simple HQET calculation yields (via the graphs in Fig. [5]) 

<rV) 
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2vr 2 



(2.15) 



(2.16) 



Solving the differential equation ( 2.15Q by standard methods, we find 



CM 



12tt 



a s (m b ) 



12/25 



(2.17) 



We now turn to the renormalization of the term C^{ii)rh? c {ji)hh. The coefficient C^/z) 
runs both because of the renormalization of the operator 

3 =m 3 Jih, (2.18) 

and because of the mixing from Oq. It obeys the renormalization group equation 

^ ° M = 73 ° M + 763 ° M ' (2,19) 

Since the current h h is conserved, the anomalous dimension 73 is given solely by the renor- 
malization of the mass m c (/x), 



73 = 37„ 



3gM 
2tt 2 



The mixing of Oq with O3 is given by contracting the charm quark fields as in Fig. 

3 



763 



2tt 



ml . 



2 ""b 



[2.20) 



(2.21) 



Because the graph has a cubic divergence, it is proportional to m 3 ; hence 763 is proportional 



to m|. 



Solving the differential equation ( [2.19|) for C 3 (/i) and setting /1 = m c , we find 



C 3 {m c ) 



m 



144vr 



m 



192tt 3 c 23a s (m c ) 



^12/25 _ 2--H/25 



where 



a s (m b ) 



(2.22) 



^2.23) 



Note that because the numerator in the expression (|2.22j ) vanishes as z — > 1, the solution 
for C*3(/i) is well behaved in the limit a s (m c ) — > 0. Expanding C^(fj) in powers of a s (m c ), 
we find 



C-Mc 



-24 lnm c -48 



a s (m c ) 



7i 



In 2 m c + . . . 



(2.24) 
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For a s (m c ) = 0.41 and rh c = 0.37 (for which a s (rrib) = 0.27 and z = 1.54), we find that the 
contribution of this operator to the total rate is of order one: 

C 3 (m c ) = IVo [1-22 - 0.21 + ...] = T vo [0.99] . (2.25) 

The resummed logs change the total decay rate by 18%. In terms of the quantities rh c (mb) 
and a s (m b ) renormalized at the scale \i = m b , we find 

T = IVo jl— 2rn c (m fe ) — %m 2 {m b ) — 18m^(m fe ) mm, 

+ 23^) A «W^ /25 - 1 ) + -} 

Note that although we include the term — 18rh 3 (m b )\nm c in this renormalization group 
improved expression, there are also uncomputed terms of the same order from the two-loop 
renormalization of the dimension six operators. 



(2.26) 
(2.27) 



III. DECAYS VIA A LEFT-HANDED CURRENT 

We now apply the same analysis to the physical situation of decays mediated by the 
left-handed current J M = &7 M (1 — 7 5 )c. Including the coupling GpV cb / y/2, the total decay 
rate is related to the forward scattering amplitude via 



Gi I V c b 1 2 



J dqe-^(B\T{J^(x),r(0)}\B) xi(^ 



The tree level rate is r = G 2 F \ V cb \ 2 m b /192Tr 3 . We now expand the time ordered product in 
an operator product expansion, as before. At tree level, and for dimension n < 6, we find 



2m 



B 



hh — %rh 2 {m b )hh + ... 



(3.2) 



where the OPE is performed at the renormalization scale \x = m b . The HQET field h(x) 
is given in Eq. ( |2.3j ) , and the lowest order matrix element of h h in Eq. ( |2.5| ) . The ellipses 
denote charm-independent dimension five operators such as h(iD) 2 h, which do not induce 
terms proportional to m™. As before, the combination rh 2 hh does not run below // = m^; 
expanding rh c {m b ) in terms of m^ ole , one finds at leading logarithmic order 



im c [m b ) 



KmP° le ) 2 ^ 24/25 



m. 



pole \ 2 



32 as K)lnrn c + -a s K) 2 ln 2 m c + ... 



(3.3) 



where z = a s (m c )/a s (m b ). Hence we reproduce the order a s correction from I^ 1 ) ( |1.3|) , and 
then extend this result to resum all logarithms of the form rh 2 a™ ln n rh c in the coefficient 
functions r^ ra ^. (This constraint on the m 2 a s \ogm c was also noted by Nir 0). 

At dimension six, four quark operators of the form hTihcTzC arise, just as in the case of 
vector decays. In principle, these operators could induce terms of order m 3 In m c in the total 
decay rate, when they mix with hh. Expanding the operator product, applying the Fierz 
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transformation, and dropping parity-odd operators which cannot contribute to the forward 
matrix element, we find the dimension six operators 



32tt 



2m b 



-\hhcfc - -hj"j 5 hcj"j 5 c - 3hT a hcT a f)c - hT a ^ 5 hcT a ^ 5 c 
2 6 



(3.4) 



Note the absence of a term proportional to hhcc, the only dimension six operator which 
can mix with m\h h. Hence, unlike the case of vector decay, there is no term in the decay 
rate proportional to w? c lnm c . In fact, since none of the dimension six operators of the form 
hThcTc and hT a ThcT a Tc mix with each other (as discussed in the previous section), there 
are no terms in the decay rate proportional to m^a™ ln n+1 m c , for any n. The absence of such 
logarithmic terms at tree level in r^ ^ is thus extended to all orders. The leading logarithms 
at order m\ are therefore simply resummed by replacing mJ? ole with m c (irib) in the expression 
forrW, 

— m c {m b f = — (mP° le ) 3 ^ 36/25 = — « olc ) 3 {l + 6a s (m b ) lnm + ... } . (3.5) 



To reproduce the m^lnm c term in r^, we must continue the OPE to include operators 
of dimension seven. Although there is a large number of such operators, we can use heavy 
quark symmetry and the classical equations to motion to reduce these to only a few which 
are relevant to the analysis. There are three classes of dimension seven operators, each of 
which can have a singlet or octet color structure and one of the four Dirac structures ( |2.12 ) . 



The first class is dimension six operators multiplied by an additional factor of m c . Because 
we are counting powers of 1 we will treat these operators as dimension seven. These 
operators do not mix with each other, just as their dimension six counterparts do not. Hence 
the only one of these operators which can mix with m\hh is mjihcc. 

The second class of operators is those in which a derivative acts on the charm quark. 
We may use the classical equation of motion Upc(x) = m c c(x) to reduce some of these to 
operators of the first class. Of those that remain, only hhafiv -iDc mixes with m 4 /i h, via the 
graphs in Fig. |4|. However, these operators can mix with each other under renormalization, 
as well as with those of the first class, via the diagrams in Figs. [5] and []. Let us consider the 
gauge in which v ■ A — 0, where the gluon does not couple to the heavy quark field. Then we 
have only the graphs which renormalize the charm quark part of the operator. These can 
mix neither Dirac nor color structures, nor can they induce operators in which a derivative 
acts on the bottom quark. The running of these operators is diagonal as well. Finally, of 
the operators of this class, only hhcfv ■ iDc mixes with rhhhcc. 

The third class is those operators in which a derivative acts on the bottom quark. Here 
the classical equation of motion v ■ iDh(x) = eliminates some operators entirely. The 
remaining ones cannot mix directly with rh 4 h h at order but they can mix with the two 
other dimension seven operators we have identified so far. In v ■ A = gauge, such mixing 
can only occur via the first two graphs in Fig. ^ Only color octet operators can mix with 
the color singlet h hc^v ■ iDc by such one-gluon exchange. The operators in this class can 
also mix among themselves. However, it is straightforward to show that while color octets 
mix with both singlets and octets, color singlets mix only with each other. Hence, in this 
class of operators, only the color octets are relevant to our calculation. Furthermore, the 
Dirac structure is sufficiently constraining that of these, only hT a \D^hcT a r )^ 1 c contributes. 
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b bb bb bb b 

FIG. 6. Additional diagrams which renormalize O3, O4 and O5. 
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There is one other source of dimension seven operators which we must consider. The 
HQET lagrangian beyond leading order is given by [18,19 



where 



£ H qet = hv ■ iDh + —O k + — C g (ji)O g (h) + . . . , (3.6) 



O k =- h(iB) 2 h and O g = | ha^G^h (3.7) 



are the leading spin and flavor symmetry violating corrections to the nib — > oo limit. The 
operators Ok and Og are treated as perturbations; because they come with explicit factors 
of 1 /mb, they can induce mixing between operators at different order in the l/m& expan- 
20| . In particular, they can mix operators of dimension six with those of dimension 



sion 

seven, via the graphs shown in Fig. [7[ We find that insertions of Og do not induce mixing 
with any of the dimension seven operators of interest, while Ok does, if the dimension six 
operator which is being renormalized is hT a hc^T a c. 

With this taxonomy in hand, we can now write down the list of operators of dimension 
six and seven which are relevant for our analysis. In the end, it is mercifully short: 

O x = hT a h cfT a c 

2 = hT a h q^T a qi , q { = u, d, s 
= ^-h{-i%T a + T a i%}h crfT a c 

04 = -—h{-\D^T a + TaVfi^h q^Taqi , q { = u, d, s (3.8) 

IVTLb 

5 = + i%}c 

ZlTLb 

n ra c (/i) -r 

Oq = till cc 

m b 

The operators O3, O4 and O^ have been constructed to be Hermitian. We also must include 
the quark bilinear m\hh. It is convenient to define it with an inverse factor of the strong 
coupling, because this will make the anomalous dimension matrix homogeneous in g 2 : 

07= i<M- hh (39) 

Factors of l/m& have also been included in O3, . . . ,0? so that the anomalous dimension 
matrix will have no explicit factors of l/m^Q To summarize, the operators 0\ through 
Oq are renormalized via the graphs in Fig. [|. In addition, the dimension seven operators 
O3, . . . ,Oq get contributions from the graphs in Fig. ^. The dimension six operators 0\ and 
O2 mix with each other via the "penguin" diagrams in Fig. §, as do O3 and O4. Finally, 0\ 
mixes with C 5 and Oq via time-ordered products with Ok, as shown in Fig. |7|. 



2 We did not include these factors in the analysis of the vector decay, because the renormalization 
group equations were already so simple. 
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FIG. 8. "Penguin" diagrams mixing 0\ and 0<i- The same diagrams mix O3 and O4. 



The operators in Eq. ( p.8|) are in fact not all independent; (D\ and O3 (and O2 and 
O4) are related via reparameterization invariance pT]]. Since f M and must appear in the 
combination 

V M = ^ + i^/2m«, - z rJ M /2m 6 , (3. 10) 

we find the restrictions 

C 1 (/x) = C 3 (/i), C 2 (/x) = C 4 (/i) . (3.11) 

Our explicit calculations confirm this result. 

We now perform the operator product expansion at the scale = m b . For the operators of 
dimension seven, the matching coefficients are generated at subleading order in the expansion 
in l/m&. The momentum of the heavy b quark is written as p^ = m^v 1 * + k^, where /c M 
is the "residual" momentum. For an on-shell h field, the classical equation of motion is 
v ■ k — H . The expression of the b quark spinor u b in terms of the heavy spinor Uh is also 
affected, becoming u b = (1 + fi,/2m b )uh . 



There are then two sources of matching onto operators of dimension seven. First, the 
lepton momentum g M may be written as g M = p% — p£ = m b v^ + k^ — The momenta k^ 
and p% lead to operators with covariant derivatives acting on the h and c fields, respectively. 
Second, the correction to the heavy quark spinors must be accounted for. We define reduced 
operator coefficients by 

Q(M) = 327 f 9 F ° Q(A«). (3-12) 

m£ 2m# 

Then performing the operator product expansion at tree level yields the nonzero terms 

2 - 1 

d(m b ) = -3 , C 3 (m b ) = -3 , C 5 (m b ) = - , C 6 (m b ) = - . (3.13) 
The Wilson coefficients evolve according the renormalization group equation 

V^piifi = 7iAG«) • (3-14) 
The anomalous dimension matrix 7^ is defined by the operator renormalization 
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^■•r (n) 



' d „ d d \ , n) 

^ + + 7m c m c - n7 cxt 



(3.15) 



where Fq} is an n-point Green function with a single insertion of the operator Cj. Using 
the known mass and wavefunction anomalous dimensions P,3,|Bi|Tl| 



7"i c 



2tt 2 



7c 



7fe 



the Feynman diagrams in Figs. 



12vr 2 ' " l 6tt 2 ' 
yield the anomalous dimension matrix 



7y 



7T^ 



/ -23/48 


1/12 








-2/9 


-1/9 \ 


1/4 


-5/16 




















-23/48 


1/12 





1/6 








1/4 


-5/16 




















4/9 


-1/9 -3/8 

















-3/2 


V o 














23/24 / 



(3.16) 



(3.17) 



The renormalization group is then used to evolve the coefficients Ci(ji) from /i = m b to 
ji = m c . The logarithmically enhanced terms m\a n s ln n+1 rh c are given by the combination 



777, f 777 1 

F(m c ) = C 7 (m c ) c[ c) (B\hh\B) 
g z (m c )m b 



^ -4A s r 



a s (m c ) 



mtC 7 (m c )-^(B\hh\B) . (3.18) 



2m 



B 



By inspection of the matching coefficients (|3.13|) and the anomalous dimension matrix ( p.!7[ ), 
we see that only the linear combination (9 5 — 40 6 mixes into 7 , and that this linear 
combination of operators does not run in the leading logarithmic approximation. Therefore, 
the solution to the renormalization group equation for CV(/i) is particularly simple. With 
the matrix element fl2.5|) , we find 



F(m c ) 



8tt 



a s (m c ) 



c 23 V 



z 



-23/25 > 



(3.19) 



When expanded and written in terms of a s (mj) and the reduced pole mass m c , the first two 
terms of the expression 

T 1 F(m c ) = -2Aml In m c - 96^^ In 2 m c + — ( * m \ \£ rh c 

ix 3 \ i ) 



- 12 



a s [m b ) 



7T 



m c In m c + 



104 / a s (m b ) 



7T 



m c In m c + . . . 



(3.20) 



reproduce the known tree level and one loop results from Eqs. ( |1.2| ) and (|1.3j), respectively. 
In addition, all of the leading logarithms have been resummed. 

As in the previous section we use a s (m c ) = 0.41 and m c = 0.37 to find 



F(m c ) = T [0.45 - 0.15 - 1.4 • 10~ s + ...]= T [0.30] 



(3.21) 



16 



The re-expanded contributions from the terms m\hh and m\hh are 

-^rhl{m b )(B\hh\B) = T [-1.10 - 0.37 - 0.12 + .. .] = T [-1.67] (3.22) 

and 

p _ 

— — m^(m 6 )(S | = T [0 + 0.92 + 0.46 + ...] = T [1.71] (3.23) 
2mB 

From these expansions, a large contribution at 0(a 2 mllnm c ) can be seen. As discussed in 
the introduction, one can expect a cancelation of this term by the a 2 s m^ In m c term, which 
has not been calculated here. 



IV. SUMMARY 

We have studied the operator product expansion for the process b — > civ, to understand 
better the origin of the "phase space" logarithms which appear in the total decay rate. Af- 
ter extracting the known tree level term, we have extended the analysis to include radiative 
corrections.^ In particular, we have used a renormalization group analysis to resum the lead- 
ing logarithms of the form m^a™ ln n m c , m^a™ +1 ln n m c and m\cx n s ln n+ m c . Unfortunately, 
these terms do not dominate, in any limit of the theory, over certain others which have 
been omitted. Hence the results of this calculation cannot be used to extract any reasonable 
estimate of the true size of the higher order corrections. 

The point of this calculation lies rather in the insight which it affords us into the origin 
of these logarithms, which even though not divergent, reflect sensitivity to physics which is 
far in the infrared with respect to the scale of the decaying b quark. We have exploited this 
separation of scales to resum to all orders a certain subset of the phase space logarithms. In 
so doing, we have explored more generally their relation to other logarithms which appear 
in the theory, such as the "hybrid" anomalous dimensions of the heavy weak current [|3],|J. 
The hybrid anomalous dimensions are also not numerically dominant, but by studying them 
one may investigate interesting questions of principle in the Heavy Quark Effective Theory. 
The analysis and resummation which we have performed here should be viewed in much the 
same spirit. 
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3 The one loop radiative corrections to the operator product expansion for nonleptonic D decays 



were calculated in Ref. [23], although no logarithms were resummed. 
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